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Abstract 

The channeling effect of low energy ions along the crystallographic axes and 
planes of Nal(Tl) crystals is discussed in the framework of corollary investigations 
on WIMP Dark Matter candidates. In fact, the modeling of this existing effect 
implies a more complex evaluation of the luminosity yield for low energy recoiling 
Na and I ions. In the present paper related phenomenological arguments are 
developed and possible implications are discussed at some extent. 
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1 Introduction 

It is known that ions (and, thus, also recoiling nuclei) move in a crystal in a differ- 
ent way than in amorphous materials. In particular, ions moving (quasi-) parallel 
to crystallographic axes or planes feel the so-called "channeling effect" and show an 
anomalous deep penetration into the lattice of the crystal [1, 2, 3]; see Fig. 1. 
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Figure 1: Simplified schema of the channeling effect in the Nal(Tl) lattice. The axial 
channeling occurs when the angle of the motion direction of an ion with the respect 
to the crystallographic axis is less than a characteristic angle, \l/ c , depicted there (see 
for details Sec. 2). Two examples for channeled and unchanneled ions are also shown 
(dashed lines). 

For example, already on 1957, a penetration of 134 Cs + ions into a Ge crystal was 
observed to a depth of about 1000 A [4], larger than that expected in the case the ions 
would cross amorphous Ge (~ 50 A). Afterwards, high intensities of H + ions at 75 
keV transmitted through thick (3000-4000 A) single-crystal gold films in the < 110 > 
directions were detected [2]. Other examples for keV range ions have been shown in 
ref. [5] where 3 keV P + ions moving into layers of 500 A of various crystals were 
studied. 

The channeling effect is also exploited in high energy Physics e.g. to extract high 
energy ions from a beam by means of bent crystals or to study diffractive Physics by 
analysing scattered ions along the beam direction (see e.g. ref. [6]). 

Recently [7] it has been pointed out the possible role which this effect can play 
in the evaluation of the detected energy of recoiling nuclei in crystals, such as the 
Nal(Tl) 2 . 

In fact, the channeling effect can occur in crystalline materials due to correlated 
collisions of ions with target atoms. In particular, the ions through the open channels 
have ranges much larger than the maximum range they would have if their motion 
would be either in other directions or in amorphous materials. Moreover, when a 
low-energy ion goes into a channel, its energy losses are mainly due to the electronic 
contributions. This implies that a channeled ion transfers its energy mainly to electrons 
rather than to the nuclei in the lattice and, thus, its quenching factor (namely the ratio 
between the detected energy in keV electron equivalent [keVee] and the kinetic energy 
of the recoiling nucleus in keV) approaches the unity. 

It is worth to note that this fact can have a role in corollary analyses in the Dark 

2 For completeness, it is worth to note that luminescent response for channeling in Nal(Tl) was 
already studied in ref. [8] for MeV-range ions. 
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Matter particle direct detection experiments, when WIMP (or WIMP-likc) candidates 
are considered. In fact, since the routine calibrations of the detectors are usually 
performed by using 7 sources (in order to avoid induced radioactivity in the materials), 
the quenching factor is a key quantity to derive the energy of the recoiling nucleus 
after an elastic scattering. Generally, for scintillation and ionization detectors this 
factor has been inferred so far by inducing tagged recoil nuclei through neutron clastic 
scatterings [9]; however, as it will be discussed in Sec. 3, the usual analysis carried out 
on similar measurements does not allow to account for the channeled events. A list of 
similar values for various nuclei in different detectors can be found e.g. in ref. [10]. In 
particular, commonly in the interpretation of the dark matter direct detection results 
in terms of WIMP (or WIMPdike) candidates the quenching factors are assumed to 
be constant values without considering e.g. their energy dependence, the properties 
of each specific used detector and the experimental uncertainties. An exception was 
in the DAMA/Nal corollary model dependent analyses for WIMP (or WIMPdike) 
candidates [10, 11, 12, 13] where at least some of the existing uncertainties on the qn a 
and qi values, measured with neutrons, were included. 

In this paper the possible impact of the channeling effect in Nal(Tl) crystals is 
discussed in a phenomcnological framework and comparisons on some of the corollary 
analyses carried out in terms of WIMP (or WIMP-likc) candidates [10, 11, 12, 13], on 
the basis of the 6.3 a C.L. DAMA/Nal model independent evidence for particle Dark 
Matter in the galactic halo 3 , are given. 



2 Luminosity, range and channeling in Nal(Tl) crys- 
tals 

The stopping power of an ion inside an amorphous material is given by the sum of two 
effects: its interaction with the nuclei (n) of the material and its interaction with the 
binding electrons (e)[23]. 

If E is the energy of the ion at any point x along the path, its stopping power can 
be written as: 

^L(E) = dEim ~ n ( E ) + dEi ° n ~e ( E ) (1) 
dx dx dx 

where dE '^^ n and Et ° n ~ e (jj?) are the nuclear and the electronic stopping powers of 
the ion, respectively. They can be evaluated - following the theory firstly developed 
in ref. [23] - by using some available packages, as the SRIM code [24]. 

We have to stress that the detectable light produced by a charged particle (either 
electron or ion) in scintillator detectors mostly arises from the energy loss in the 
electronic interactions. Thus, the differential luminosities in scintillators, and 

' ' dx 



3 We remind that various possibilities for some of the many possible astrophysical, nuclear and 
particle Physics scenarios have have been analysed by DAM A itself both for some WIMP/WIMP-like 
candidates and for light bosons [10, 11, 12, 13, 14], while other corollary analyses are also available in 
literature, such as e.g. refs. [15, 16, 17, 18, 19, 20, 21, 22, 7]. Many other scenarios can be considered 
as well. 
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dL J° n for electrons and ions, respectively, can be written as: 
dL e dE e _ e 

— — = a — tor electrons 

dx dx 

dLion dEi on — e . dEi on / / £ / \ /rt\ 

— - — = a (E) = a x — - — x q {Ej tor ions (recoils) (2) 

ax ax ax 

where a is a proportionality constant, dE ^~ s is the stopping power of an electron 

in the material and q'{E) = ( dE '^~° ( E )) I ( ^g ( E ) + dEi £" ( E )) is dchned as 
"differential quenching factor" . 

The total detected luminosities, L — J path ^;dx, for electrons and ions can be 
written in the form: 

L e — a [ e dx = aE e for electrons 

Jpath dx 

L ion = a j q'(E) x dE ™ n dx = ax qi(E ion ) x E ion for ions (recoils) (3) 

Jpath dx 

in addition, the range of an ion is: 



R - (E) ^ [ diEldx - «W x £ m5jdx - ^ x r ^^- 



(4) 



In eq. (3) and (4) qi (E) = ± J E q'(E')dE' is the "light quenching factor" and q 2 (E) = 

fo dE^l/dx l Jo dE~7fdi is the " ran S e quenching factor". In the energy region of 
interest for the dark matter detection q\(E) ~ q2(E) within 10-20%. In particular, 
the values of the quenching factors for recoiling nuclei in detectors made of amorphous 
materials are well below the unity in the keV energy region. 

The situation changes when the detector is either a crystal or a multi-crystalline 
material (the size of a single crystal has to be larger than few thousands of A). In this 
case the luminosity depends on whether the recoiling nucleus is (quasi-) parallel to the 
crystallographic axes or planes or not. In the first case, since the energy losses of the 
ion are mainly due to the electronic contributions, the penetration (and the range) of 
the ion becomes much larger, of the order of R e , and the quenching factor approaches 
the unity. 

The theory of ion channeling in crystals has been developed e.g. in ref. [25, 3]. In 
particular, this theory deals with channeling of low energy, high mass ions as a separate 
case from high energy, low mass ions. Here, we only remind that the condition for a 
low energy ion and a recoiling nucleus to be axially channeled along a certain string 
of atoms in the lattice is linked to a critical angle, ^> c (see Fig. 1); for details refer to 
[25, 3]. When the ion (recoiling nucleus) has a moving direction with an angle ^ with 
respect to this string lower than ^> c , it is axially channeled. 

The critical angles for axial channeling is given by [25, 3]: 
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where C 2 ~ 3 is the Lindhard's constant and d is the inter-atomic spacing in the 
crystal along the channeling direction. The Thomas- Fermi radius, citf, can be written 
as [25, 3]: 

0.8853a 

a TF = JTT (6) 

where Z\ and are the atomic numbers of the projectile (recoiling nucleus) and 
target atoms, respectively; a = 0.529A is the Bohr radius. 

The characteristic angle ^ is defined as a function of the ion (recoiling nucleus) 
energy, E, by: 



/2ZiZ 2 e 2 

where e is the electron's charge. The condition f < $ c for axial channeling is valid for 
*i > #i lim = that is for E < E Hrn = 2Z ^ 2 e 2 d ^ gi Hence? typical values for 

Nal(Tl) crystal assure that for recoil's energies less than 170 keV the quoted formulas 
hold. For completeness, we just remind that it has also been suggested that the critical 
angles may slightly depend on the temperature [26]. Moreover, the critical angles at 
low energy have a weaker dependence on ion energy than those at higher energy In 
fact, at higher energy, the critical angle is ~ [25, 3]. 

In the case of planar channeling for low energy ions the critical angle can be written 
as [27] (see also ref. [25]): 

Z!Z 2 e^ 1/3 



e pi= a TFVNd P ^^^J , (8) 

where N is the atomic number density and d p is the inter-plane spacing. The depen- 
dence of 9pi on the energy is weaker than that at higher energy, where it can be written 
as [6, 28]: 



= O.TF \J Ndp ^ 



Ecltf J 



Taking into account the critical angles for axial and planar channeling in Nal(Tl), 
we have calculated by Monte Carlo method the solid angle interested by both axial 
and planar channeling in Nal(Tl) crystals as a function of the energy of the recoiling 
nuclei, Er; see Fig. 2. Moreover, just the lower index crystallographic axes and planes 
have been considered, for the axial channeling: < 100 >, < 110 >, < 111 > and for 
the planar channeling: {100}, {110}, {111}. 

In this way, the estimated light response of a Nal(Tl) crystal scintillator to Sodium 
and Iodine recoils at given energy has been studied taking into account the channeling 
effect in the considered modeling. For a given nucleus A with recoil energy Er the 
response of a Nal(Tl) crystal scintillator can be written as (Edet\En), where 

Edet is the detected energy. By the definition: J °° ^Eat (Edet\En)dEd e t = 1- In 
most cases of the Dark Matter direct detection field - that is without including the 
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Figure 2: Fraction of solid angle interested by both axial and planar channeling in 
Nal(Tl) crystals as a function of the energy of the recoiling nuclei, calculated according 
to the text. In these calculations just the lower index crystallographic axes and planes 
have been considered: for the axial channeling: < 100 >, < 110 >, < 111 > and for 
the planar channeling: {100}, {110}, {111}. 

channeling effect - the light response is assumed equal to a Dirac delta function: 
^ (Ed et \Efi) = 8(Edet — QaEr), where qa is the value (assumed constant) of the 
quenching factor of the unchanneled A nucleus recoils. 

The evaluation of , when accounting for channeling effect, has been realised by 
means of a Monte Carlo code; the path of a recoiling nucleus, at a given recoil energy 
Er, has been calculated under the following reasonable and cautious assumptions: 

i ) isotropic distribution of the recoils; 

ii ) in the case the recoil would enter in a channel, a de-channeling can occur due to 
some interactions with impurities in the lattice, as Tl luminescent dopant centers. 
The probability density of such a process is assumed to be: p(x) = je~ x / x , with 
A = 1200 A, that is the average distance among the Tl centers; 

iii ) the energy losses by the recoil nuclei in a channel just depend on the electronic 
stopping power (see eq. (1)); 

iv ) the energy losses by the recoil nuclei in a channel arc converted into scintillation 
light with a quenching factor ~ 1; 

v ) if a recoil is de-channeled due to a nuclear interaction, it can either re-enter into 
another channel or not. The differential distribution of the nuclear interaction 
is assumed to be isotropic; 

vi ) in case the recoil would not enter into a channel, it is cautiously assumed 
that it stops and its released energy is converted into scintillation light with the 
quenching factor of unchanneled events. In this case the straggling is considered 
as evaluated by the SRIM code [24]. 



G 



0.25 



0.2 



0.15 



0.1 



0.05 



0.2 
0.175 
0.15 
0.125 
0.1 
0.075 
0.05 
0.025 




a) Iodine 
E R =4 keV 



E„ et (keVee) 

c) Sodium 
E R =4 keV 



0.014 
0.012 
0.01 
0.008 
0.006 
0.004 
0.002 


x 10 ' 

0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 




b) Iodine 
E R =40 keV 



10 
E. 



20 30 



40 



^(keVee) 

d) Sodium 
E R =40 keV 



E det (keVee) 




20 30 
(keVee) 



Figure 3: Examples of light responses in terms of keVee, ^E^t (^det\En), for Iodine 
recoils of 4 keV (a) and of 40 keV (6) and for Sodium recoils of 4 keV (c) and of 40 
keV (d) in the modeling given in the text. In this calculation the quenching factors 
for Sodium and Iodine recoils in amorphous or out of channel Nal(Tl) are assumed at 
the mean values given in ref. [29]. Just to emphasize the effect of the channeling, the 
broadening due to the energy resolution of the detector has not been included here. 
The peaks corresponding to fully channeled events (q ~ 1) and to fully quenched events 
(broadened by the straggling) are well evident; in the middle events, which have been 
de-channeled at least once, are also visible. It is possible to note that e.g. in the case 
of Iodine recoils the fully channeled events are about 25% at 4 keV; this percentage 
becomes smaller, about 1% at 40 keV. 



In Fig. 3 few examples of light responses in terms of keVee for Iodine recoils of 4 keV 
(a) and of 40 keV (b) and for Sodium recoils of 4 keV (c) and of 40 keV (d) are given. 
In these calculations the quenching factors for Sodium and Iodine recoils in amorphous 
or out of channel Nal(Tl) are assumed at the mean values given in ref. [29]. Just to 
emphasise the effect of the channeling, the broadening due to the energy resolution of 
the detector has not been included in this figure. The peaks corresponding to fully 
channeled events (q ~ 1) and to fully quenched events (broadened by the straggling) 
are well evident; in the middle events, which have been de-channeled at least once, are 
also visible. 

Finally, in Fig. 3 it is possible to note that the number of fully channeled (q ~ 1) 
events strongly decreases when increasing the recoil energy: they are ~ 25% at 4 keV 
and ~ 1% at 40 keV for Iodine recoils and ~ 18% at 4 keV and ~ 0.3% at 40 keV for 
Sodium recoils. These behaviours are depicted in Fig. 4. 
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E R (keV) 

Figure 4: Fraction of events with quenching factor ~ 1, that is fully channeled events, 
as a function of the energy of the recoiling nuclei in Nal(Tl) crystals according to the 
modeling described in the text. 

3 Some comments 

Let us now analyse the phenomenologies connected both with the data on nuclear 
recoils induced by neutron scatterings and with the WIMP (or WIMP-like) direct 
detection in the light of the presence of the channeling effect. 

In particular, Fig. 5 shows some examples of neutron calibrations of Nal(Tl) de- 
tectors at relatively low recoil energy. There the energy responses of the used Nal(Tl) 
detectors to Sodium recoils of 10 keV [30] and of 50 keV [31] are reported as solid 
histograms; the peaks corresponding to the quenched events arc well clear. The super- 
imposed continuous curves have been calculated as those of Fig. 3, obviously broaden- 
ing them by the energy resolution of the corresponding detector. The fully channeled 
peaks (q ~ 1), which in these cases can contain only the 6% and 0.15% of the events 
respectively (see Fig. 4), are smeared out by the energy resolution and only contribute 
to the higher energy tails in the energy spectra. 

Thus, the simple analysis of Fig. 5 shows that the neutron data can contain channeled 
events, but - owing to the low-statistics of these measurements, to the small effect 
looked for and to the energy resolution - they cannot easily be identified. Moreover, 
as already shown, the channeling effect becomes less important at increasing energy 
and gives more suppressed contributions in the neutron scattering data. For Iodine 
recoils the situation is even worse. Therefore, there is no hope to single out the 
channeling effect in the already-collected neutron data. 

On the other hand, the accounting of the channeling effect can give a significant 
impact in the sensitivities of the Dark Matter direct detection methods when WIMP 
(or WIMP-like) candidates are considered. In particular in par. 4 and 5 we will show 
that lower cross sections are explorable in given models for WIMP and WIMP-like 
candidates by crystal scintillators, such as Nal(Tl) (up to more than a factor 10 in 
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Figure 5: Examples of neutron calibrations of Nal(Tl) detectors at low recoil energy. 
In particular, the energy responses of Nal(Tl) detectors to Sodium recoils of 10 keV 
(left panel) [30] and of 50 keV (right panel) [31] are shown; the peaks corresponding 
to the quenched events are well clear. The superimposed continuous curves have been 
calculated as those of Fig. 3, obviously broadening them by the energy resolution of 
the corresponding detector. The fully channeled peaks (q ~ 1), which contain in these 
cases only the 6% and 0.15% of the events respectively, are smeared out by the energy 
resolution and can just contribute to the higher energy tails in the energy spectra. 



some mass range). Moreover: 



i ) similar situation holds for purely ionization detectors, as Ge (HD-Moscow - like 
experiments); 

ii ) the loss of sensitivity occurs when pulse shape discrimination is used in crystal 
scintillators (KIMS); in fact, the channeled events (q ~ 1) are probably lost; 

iii ) no enhancement can be present in liquid noble gas experiments (DAMA/LXe, 
WARP, XENON, ...); 

iv ) no enhancement is possible for bolometer experiments; on the contrary some 
loss of sensitivity is expected since events (those with q ion ~ 1) are lost by 
applying some discrimination procedures, based on q ion « 1. 



4 Application to the WIMP-nucleus elastic scatter- 
ing 

Let us now consider the case of WIMP (or WIMP-like) elastic scattering on target 
nuclei. In particular, the expected differential counting rate of recoils induced by 
WIMP-nucleus elastic scatterings has to be evaluated in given astrophysical, nuclear 
and particle physics scenarios, also requiring assumptions on all the parameters in- 
volved in the calculations and the proper consideration of the related uncertainties 
(for some discussions see e.g. [10, 11, 14, 12, 13]). Hence, the proper accounting for 
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the channeling effects must be considered as an additional uncertainties in the eval- 
uation of the expected differential counting rate. The usual hypothesis that just one 
component of the dark halo can produce clastic scatterings on nuclei will be assumed 
here. In addition, the presence of the existing Migdal effect and the possible SagDEG 
contribution - we discussed in refs. [13] and [12] respectively - will be not included 
here for simplicity. Thus, for every target specie A, the expected distribution of the 
detected energy can be written as a convolution between the light response function, 
dE^ t i defined in the previous section, and the differential distribution produced in the 
WIMP-nucleus elastic scattering: 



d^A i n f dN A ,dR A 



(oh) 

, r (E det ) = I ^{E det \E R )^{E R )dE R . (10) 
dhdet J dhdet dh R 

The differential energy distribution of recoils, as function of the recoil energy E R , is: 

§±{E R )=N T ™- P"" -^(v,E R )vf(v)dv. (11) 

In this formula: i) is the number of target nuclei of A specie; ii) pw = £po, where 
p is the local halo density and £ < 1 is the fractional amount of local WIMP density; 
hi) mw is the WIMP mass; iv) f(v) is the WIMP velocity (v) distribution in the Earth 
frame: v) v m ; Ln = A / ™ A : i ER (via and raw a are the nucleus mass and the reduced mass 

V WA 

of the WIMP-nucleus system, respectively) ; vi) v max is the maximal WIMP velocity 
in the halo evaluated in the Earth frame; vii) -M-(v,E R ) = (-M-' ^ ' da 



'It:,: - \dE R ) SI 1 \dE RJsD 



with ( -jj^ j spin independent (SI) contribution and ( \ spin dependent (SD) 
contribution. 

Finally, the expected differential counting rate as a function of the detected energy, 
Edet, for a real multiple-nuclei detector (as e.g. the Nal(Tl)) when taking into account 
the channeling effect can easily be derived by summing the eq. (10) over the nuclei 
species and accounting for the detector energy resolution: 

^l(E det )= G(E deU E') ]T ^(E')dE'. (12) 

The G(E det , E') kernel generally has a gaussian behaviour. 

Few examples of shapes of expected energy distributions with and without account- 
ing for the channeling effect, calculated in the modeling given above, are shown in Fig. 
6. For this template purpose - accounting also for the experimental features of the 
detectors [32, 33, 10, 11] - we have just adopted the following additional assumptions 
(among all the many possibilities): i) WIMP mass of mw — 20 GeV; ii) WIMP with 
dominant Spin Independent coupling and with nuclear cross sections cx A 2 ; iii) point- 
like SI cross section a si = 10~ 6 pb; iv) an halo model NFW (identifier A5 in ref. [34], 
local velocity vo = 220 km/s and halo density at the maximum value 0.74 GeV cm~ 3 
[34]; v) form factors and quenching factors of 23 Na and 127 I as in case A of ref. [10]. 
These pictures point out the enhancement of the sensitivity due to the channeling 
effect according to the given modeling. 
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Figure 6: An example of the shapes of the expected energy distributions in Nal(Tl) 
from Sodium and Iodine recoils induced by WIMP interactions with (continuos line) 
and without (dashed line) including the channeling effect in the crystal for the scenario 
given in the text. Left panel: behaviour of the unmodulated part of the expected signal, 
So. Right panel: behaviour of the modulated part of the expected signal, S m - The 
vertical lines indicate the energy threshold of the DAMA/Nal experiment. 



5 Examples of the possible impact on some corollary 
quests from the DAMA/Nal data 

The accounting of the channeling effect in corollary quests for WIMPs as Dark Matter 
candidate particles can be investigated by exploiting the expected energy distribution, 
derived above, to some of the previous analyses on the DAMA/Nal annual modulation 
data in terms of WIMP-nucleus elastic scattering. For this purpose, the same scaling 
laws and astrophysical, nuclear and particles physics frameworks of refs. [10, 11] are 
adopted. In addition, as already mentioned, for simplicity just to point out the impact 
of the channeling effect, the possible SagDEG contribution to the galactic halo and 
the presence of the existing Migdal effect - whose effects we discussed in refs. [12] and 
[13], respectively - will not be included here. 

The results for each kind of interaction are presented in terms of allowed vol- 
umes/regions, obtained as superposition of the configurations corresponding to likeli- 
hood function values distant more than 4<r from the null hypothesis (absence of mod- 
ulation) in each one of the several (but still a very limited number) of the considered 
model frameworks. This allows us to account - at some extent - for at least some of the 
existing theoretical and experimental uncertainties (see e.g. in ref. [10, 11, 14, 12, 13] 
and in the related astrophysics, nuclear and particle physics literature). Here only the 
low mass volumes/regions, where the channeling effect is dominant, are depicted. 

Since the 23 Na and 127 I are fully sensitive both to SI and to SD interactions, the 
most general case is defined in a four-dimensional space (myy, £&sij ti&SD, 0), where: 
i) a si is the point-like SI WIMP-nucleon cross section and ctsd is the point-like SD 
WIMP-nucleon cross section, according to the definitions and scaling laws considered 
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in ref. [10]; ii) tgO is the ratio between the effective coupling strengths to neutron and 
proton for the SD couplings (9 can vary between and ir) [10]. 




m w (GeV) 

Figure 7: An example of the effect of the channeling, modelled as in the text, on a 
DAMA/Nal allowed region for purely SI coupling WIMPs in the given scenario (see 
text). The region encloses configurations corresponding to likelihood function values 
distant more than 4<r from the null hypothesis (absence of modulation) . This example 
has been evaluated according to the assumptions given in the text. In particular, an 
halo model Evans' logarithmic with R c = 5 kpc (identifier Al in ref. [34]) has been 
considered for a vq value of 170 km/s and halo density at the corresponding maximum 
value [34]; the form factors parameters and the quenching factors of 23 Na and 127 I are 
as in case A of ref. [10]. The solid (dashed) curves delimitate the allowed regions when 
the channeling effect is (not) included. For simplicity just to point out the impact of 
the channeling effect, the possible SagDEG contribution to the galactic halo and the 
presence of the existing Migdal effect - whose effects we discussed in refs. [12] and 
[13], respectively - are not included here. Moreover, the same considerations reported 
in ref. [10] still hold. 

Preliminarily, here to offer an example of the impact of accounting for the channel- 
ing effect as given in the text, Fig. 7 shows a comparison for allowed slices correspond- 
ing to purely SI coupled WIMPs in some particular given scenario. This example has 
been evaluated for an halo model Evans' logarithmic with R c = 5 kpc (identifier Al in 
ref. [34]) for a vo value of 170 km/s and halo density at the corresponding maximum 
value [34]; the form factors parameters and the quenching factors of 23 Na and 127 I are 
as in case A of ref. [10]. The solid (dashed) curves delimitate the allowed regions in 
the given scenario when the channeling effect is (not) included. As it can be seen, 
for WIMP masses in the few-20 GeV region the allowed SI region when including the 
channeling effect is lower than one order of magnitude in cross section. 

The subcase of purely SI coupled WIMPs for the scenarios of ref. [10, 11] is shown 
in Fig. 8, while in Fig. 9 four slices of the 3-dimensional allowed volume (mw, £&sd, 
9) for the purely SD case are given as example; the low mass region of interest for the 
effect is just focused here. 

Finally, in the general case of mixed SI&SD coupling one gets, as mentioned above, 
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m w (GeV) 

Figure 8: Region allowed in the (£crsi, mw) plane in the same model frameworks of ref. 
[10,11] for pure SI coupling; just the low mass part of interest for the channeling effect 
is focused here. The dotted region is obtained in absence of channeling effect [10,11], 
while the dashed one is obtained when accounting for it as described in the text. The 
dark line marks the overal external contour. It is worth to note that the inclusion 
of other contributions and/or of other uncertainties on parameters and models, such 
as e.g. the possible SagDEG contribution [12] and the Migdal effect [13] or more 
favourable form factors, different scaling laws, etc., would further extend the region 
and increases the sets of the best fit values. For completeness and more see also [10-14]. 
The same considerations reported in the caption of Fig. 7 hold. 
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Figure 9: Examples of slices of the 3-dimcnsional allowed volume {^(Jsd, mw, 0) in the 
same model frameworks of ref. [10,11] for pure SD coupling; just the low mass part 
of interest for the channeling effect is focused here. Analogous comments and remarks 
as those in the captions of Figs. 7 and 8 hold. 
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a 4-dimensional allowed volume (£o\sj, £,&sd, mw, 0). New allowed volume at the given 
C.L. is present in the GeV region when accounting for the channeling effect. Fig. 10 
shows few slices of such a volume as examples. 



= = 7t/4 = 71/2 = 2.435 




^a SD (pb) 

Figure 10: Examples of slices of the 4-dimensional allowed volume (£<jsi, £,(tsd, mw, &) 
in the model frameworks considered in ref. [10,11] for mixed SI&SD coupling; just the 
low mass part of interest for the channeling effect is focused here. Analogous comments 
and remarks as those in the captions of Figs. 7 and 8 hold. 

Note that general comments, extensions, etc. already discussed in ref. [10, 11, 14, 
12, 13] still hold. 

6 Conclusions 

In this paper the channeling effect of recoiling nuclei induced by WIMP and WIMP- 
like elastic scatterings in Nal(Tl) crystals has been discussed. Its possible effect in a 
reasonably cautious modeling has been presented as applied to some given simplified 
scenarios in corollary quests for the candidate particle for the DAMA/Nal model in- 
dependent evidence. This further shows the role of the existing uncertainties and of 
the correct description and modeling of all the involved processes as well as their pos- 
sible impact in the investigation of the candidate particle. Some of them have already 
been addressed at some extent, such as the halo modeling [34, 10, 11], the possible 
presence of non-thermalized components in the halo (e.g. caustics [35] or SagDEG 
[12] contributions), the accounting for the electromagnetic contribution to the WIMP 
(or WIMP-like) expected energy distribution [13], candidates other than WIMPs (e.g. 
[14] and in literature), etc.. 
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Obviously, many other arguments can be addressed as well both on DM candidate 
particles and on astrophysical, nuclear and particle physics aspects; for more see [10, 
11, 14, 12, 13] and in literature. In particular, we remind that different astrophysical, 
nuclear and particle Physics scenarios as well as the experimental and theoretical 
associated uncertainties leave very large space also e.g. for significantly lower cross 
sections and larger masses. 
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